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ABSTRACT 

In a relativistic model of a neutron star, the star's mass is less than the mass of the individual 
component baryons. This is due to the fact that the star's negative binding energy makes a contribution 
to the star's total energy and its mass. A consequence of this relativistic mass deficit is that a neutron 
star that is accreting matter increases its mass at a rate which is slower than the mass of a baryon 
times the rate that baryons are accreted. This difference in the rate of change of the masses has a simple 
relation with the star's gravitational redshift. We show that this effect has the potential to be observed in 
binaries where the mass transfer is driven by angular momentum losses from the gravitational radiation 
emitted by the binary motion, if the physics of the donor star is well-enough understood. 

Subject headings: stars: neutron — relativity — gravitational radiation — binaries — accretion 
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1. INTRODUCTION 

In the general theory of relativity, all forms of stress 
and energy act as a source for the gravitational field. As 
a result, the mass of a compact star is the total energy of 
the star divided by c^. As the total energy of a star in- 
cludes both the rest-mass energy and the negative binding 
energy, the mass of the star is smaller than the sum of the 
masses of its component particles. This mass deficit effect 
is roughly proportional to the star's compactness, so it is 
most important for neutron stars. For neutron stars near 
the upper mass limit, the mass deficit can be as large as 
25% of the star's mass (see for example. Cook, Shapiro & 
Teukolsky (1994)), depending on the equation of state. 

A consequence of the mass deficit is that when a neu- 
tron star accretes matter, its mass changes at a rate that 
differs from the rate that baryons are accreted times the 
mass of a baryon. This mass deficit effect is then poten- 
tially observable through changes in the orbital period of 
the binary caused by accretion, if the mass accretion rate 
is known. An important class of binaries where this effect 
may be observable are binaries in which the mass transfer 
is driven by angular momentum losses due to the emission 
of gravitational radiation caused by the orbital motion. 
In the standard scenario for gravitational-radiation-driven 
mass transfer (Kraft et al. 1962; Faulkner 1971), the mass 
deficit is neglected since the primaries are white dwarf 
stars. However, in binaries where the primary is a neu- 
tron star, the effect is in principle observable, and could 
lead to constraints on the neutron star mass and equa- 
tion of state. In this paper we will address this issue and 
discuss the conditions where the relativistic mass deficit's 
effects could be detected. The effect described requires 
very precise timing of the binary's orbital period in order 
to be detected, so we will focus our attention on the class 
of binaries containing an accreting millisecond X-ray pul- 
sar. At present three binaries of this class have measured 
mass functions: SAX J1808. 4-3658 (Wijnands & van der 
Klis 1998; Chakrabarty & Morgan 1998), XTE J1751-305 
(Markwardt et al. 2002) and XTE J0929-314 (Galloway 



et al. 2002). In addition two more pulsars, XTE J1807-294 
(Markwardt et al 2003) and XTE J1814-338 (Strohmayer 
et al 2003) have recently been discovered, although their 
mass functions are not yet measured. Precise timing of the 
X-ray pulse arrival times allows the detection of these bi- 
nary systems' orbital periods through the orbital Doppler 
effect. It is expected (Chakrabarty & Morgan 1998) that 
future observations will reveal an orbital period derivative. 

The effect of changing a neutron star's binding energy 
in a binary system has been examined by other authors in 
other situations. In the case of a binary system without 
mass-transfer, Spyrou & Stergioulas (2001) considered the 
changes in a neutron star's binding energy due to mag- 
netic dipole spin-down and its effect on the orbital pe- 
riod. The magnitude of the effect examined by Spyrou & 
Stergioulas (2001) is quite small, but may be measurable. 
In the present paper we allow for mass transfer, which 
produces an effect that is much larger. While accretion- 
induced changes in a neutron star's characteristics have 
been considered by Spyrou (1988), his model was for a sim- 
ple Newtonian stellar model and did not consider changes 
in the orbital period. 

In section 2 we review the concept of the binding energy 
of a star and discuss its value and its derivatives for rel- 
ativistic models of neutron stars. In section 3 we review 
the scenario of gravitational-radiation-driven mass trans- 
fer and include the mass deficit in the derivation. In sec- 
tion 4 we review the different theories of how the compact 
binaries of interest may have evolved and discuss which 
theories would allow the desired measurements of the rel- 
ativistic effect. Finally in section 5 we discuss the main 
theoretical uncertainties in this method. 

2. THE BINDING ENERGY OF COMPACT OBJECTS 

If we define N to be the total number of baryons in the 
star, vriB the mass of a baryon and e the binding energy 
of the star, then the total mass, Af , of the star is 

M^jubN + s/c^. (1) 
The binding energy is the sum of the star's gravitational 
potential energy, the internal energy and the rotational 
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kinetic energy. If the star were to be disassembled into 

its component baryons, the mass of the baryons would be 
Mq = uibN. For a gravitationally bound system, the 
binding energy is negative, so M is always smaller than 
Mq. Using Newtonian arguments, the leading order con- 
tribution to the binding energy is £ —GM^/R, so that 
the mass is given by 

M = Mo(l-.o^) (2) 

where Eq is a constant of order unity. Equation (2) shows 
that the relativistic correction to the mass of a star scales 
as the compactness GM/{Rc^). Hence this correction is 
negligible for main sequence stars and white dwarfs. How- 
ever, for neutron stars, the compactness can be as large as 
4/9. The contribution of the binding energy to the mass 
can be significant for compact stars. 

It is useful to briefly review the concept of mass and 
binding energy by recalling their definitions for a spherical 
star in general relativity (see box 23.1 of Misner, Thorne 
& Wheeler (1973) for an elementary discussion). A zero 
temperature equation of state (EOS) can be written as 
equations for P{n) and e{n) where n is the number density 
of baryons, P is the pressure and e is the energy density. 
The energy density includes both the matter's rest-mass 
energy density and its internal energy density. Once the 
EOS is specified, the star's total mass (using Schwarzschild 
coordinates) is 



/ 

Jo 



M = A-K i er^dr. 



(3) 



This is the mass which is measured by Kepler's law 
when a satellite orbits the star. For this reason, M 
is often called the "gravitational mass" . The baryon 
mass, Mq, of the star is given by the volume integral 
of the baryon number density times ms- Since the 
volume element in Schwarzschild coordinates is dV = 
47rr^ (1 — 2m(r)/r)~^^'^ dr, the baryon mass is 

-1/2 2, 

' r dr. 



Mq = uibN = A-ktub I n(l-2m(r)/r) 
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(4) 



Even in the case of an incompressible fluid, there is a dif- 
ference between M and Mq due to the binding energy of 
the star. These definitions can be extended to rotating 
stars (see for example, Friedman et al. (1986)). The result 
is a family of equilibrium stellar models with gravitational 
masses that depend on the star's baryon number N and 
its angular momentum J (Bardeen 1973). 

Consider the accretion of baryons with specific angular 
momentum ^ at the rate N onto a star. Since the grav- 
itational mass is a function of the stars' baryon number 
and angular momentum, the rate of increase of the gravi- 
tational mass is 

'(dM\ JdM\ 

The partial derivatives of the gravitational mass are given 
by (Bardeen 1973) 

('-^\ A, (6) 



(7) 



M = N 



(5) 



where is the star's angular velocity, and 



where $ is the dimensionless chemical potential (or injec- 
tion energy) for bringing one baryon from infinity to the 
pole of the rotating star. Using the metric of an axisym- 
metric rotating star, given by Friedman et al. (1986) 

= -e^'^dt^ + e^'l' {d(j) - ujdtf + e^i' [dr'^ + r'^dO^) , (8) 
the chemical potential is given by 

* = (Op, (9) 

where the siibscript "p" denotes that the function is eval- 
uated at the star's spin pole. In the limit of a nonrotating 
star. 



lim $ 

n->o 



1 - 



2GM 



(10) 



Typical values for the chemical potential for rotating neu- 
tron stars can be found in the paper by Friedman et al. 
(1986), however they define a quantity (3 that is related to 
the chemical potential by $ = y/]3. The chemical potential 
can be thought of as the specific energy required to bring 
a particle from infinity to the star. Since a particle at rest 
at infinity has larger energy than a particle at rest at the 
surface of the star (due to gravitational redshift), energy 
is liberated during accretion. The quantity 1 — $ is the 
fraction of the energy liberated during accretion^ 
If we define a dimensionless angular velocity, Q, 



n = n 



R^ 
GM 



and a dimensionless specific angular momentum, £, 

c£ 

i = 



MG' 



(11) 



(12) 



the rate of change of the star's gravitational mass due to 
accretion is 

M = tubN (^^ + m ' j . (13) 

In general, not all of a particle's angular momentum will 
be available to be added to the star's angular momentum. 
External torques, such as the star's magnetic field (Ghosh 
& Lamb 1991) or gravitational radiation (Wagoner 1984; 
Andersson et al. 1999) from nonaxisymmetric perturba- 
tions may remove the angular momentum. It is useful 
then to introduce a torque parameter r} defined so that if 
the baryon is in a circular orbit at radius r around the star 
when it is accreted, the specific angular momentum added 
to the star is 



(14) 



where £k is the dimensionless specific angular momentum 
of a particle in a circular Keplerian orbit. In the limit of 
a nonrotating star. 



lim t 



re 



M 



K 



GMV rc^lG-ZM' 



(15) 



For a rotating star Ik is given by the quantity cL/MG, 
where L is the specific angular momentum defined in equa- 
tion (7) of Cook, Shapiro & Teukolsky (1994). Since the 
rate that the star's mass increases will be close to uib 
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Fig. 1. — Plots of the pajrameter 7 versus neutron star mass for different models of nonrotating neutron stars. 



multiplied by the rate that baryons are accreted, we will 
introduce a parameter 7, defined by 



7=1- 



M 



3/2 



(16) 



which is roughly proportional to the star's compactness. 
The parameter 7 corresponds to the fractional rate of 
change of the relativistic mass deficit. For a slowly ro- 
tating star, 



lim 7 : 
a->o 



1 - 



2GM 



(17) 



Recall that for a non-rotating star, the gravitational red- 
shift, z, is defined by 

z = . ^ - 1, (18) 



2GM 
Rc^ 



SO in the limit of spherical symmetry, the parameter 7 

is related to z by 7 = z/{z + 1). For rotating stars, 
this relationship can be generalised by defining Zp clS 
the gravitational redshift from the star's pole. The pa- 
rameter 7 is then related to the polar redshift by 7 = 

Zp/{zp + 1) - v^K^ {GM/mf^. 

In Figure 1 plots of 7 versus gravitational mass for non- 
rotating neutron stars are shown for a few representative 
equations of state. We have chose three traditional neu- 
tron star equations of state (labeled C, L and APR) and 
two simplified quark star (labeled B60 and B80) equations 
of state. Prom the Arnett & Bowers (1977) catalogue, we 
have chosen equation of state C (Bethe & Johnson 1974) 
which is moderately stiff and L (Pandharipande & Smith 
1975) which is very stiff in order to show a typical range 
of stiffness for equations of state. Equation of state APR 
(Akmal et al. 1998) makes use of data from modern nu- 
cleon scattering experiments and also includes first order 

^ The code is available at: http://www.gravity.phys.UMm.edu/rns/ 



special relativistic corrections. The quark star equations 
of state are constructed in an approximation described in 
Glendenning (2000) where the quarks' masses are set to 
zero and the MIT bag model is used. The number after 
the letter "B" in the quark equations of state refer to the 
value of the bag constant, B, in units of McV/fm"^. 

The effects of rotation are shown in Figure 2, where 
values of 7 for the APR equation of state are plotted for 
different spin frequencies Vg and torque values. The rela- 
tivistic parameter 7 was computed by evaluating equation 
(16) after computing the value of the rotating star metric 
given in equation (8) using an accurate 2D code^ written 
by Stergioulas and described in Stergioulas & Friedman 
(1995). This code is, in turn, based on the methods de- 
scribed by Komatsu, Eriguchi & Hachisu (1989) and Cook, 
Shapiro & Teukolsky (1994). In the case where the ac- 
creted matter does not add any angular momentum to the 
star (77 = 0), it can be seen in Figure 2 that the func- 
tion 7 is almost independent of the star's spin frequency. 
The plots labeled rj — 1 correspond to stars where the 
magnetic field's effect on the accretion disk is neglected 
so that the orbits extend down to the marginally stable 
circular orbit and the accreted matter adds all of its an- 
gular momentum to the star. In this case, the effect of 
the star's spin is to decrease 7 by a large fraction over 
the nonrotating case. Since the value of has very little 
angular velocity dependence, it can be seen that the ef- 
fect of the star's angular velocity is to approximately shift 
the curve downwards by an amount linear in the angu- 
lar velocity. If the torque parameter were negative, the 
curves would be shifted vertically upwards by an amount 
linear in the angular velocity. The magnetic disk models 
of Rappaport, Fregeau, & Spruit (2004) allow material to 
be accreted without being propellered away from the star 
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Fig. 2. — Plots of the parameter 7 versus neutron star mass for equation of state APR, for different stellar spin frequencies {vg) and torque 
parameters (r;). 
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while decreasing the star's angular momentum. In their 

models, the torque parameter can be negative and have 
a magnitude which is of order unity. The curves labelled 
T] = —1 correspond to this scenario and have an increased 
value of 7. 

Simultaneous measurements of a neutron star's mass 

and the parameter 7 have the potential to constrain the 
neutron star equation of state. Since 7 is defined through 
the difference of M and tubN, accurate measurements of 
these derivatives would be desirable. Clearly these are not 
trivial measurements. In the highly idealized case of mass 
transfer driven by gravitational radiation, there is the po- 
tential to make this measurement if enough is known about 
the companion star. In the next section we review the 
gravitational-radiation-driven mass-transfer scenario, and 
discuss how 7 could be measured in an ideal situation. 
We will also discuss the uncertainties in the physics of the 
scenario that will likely limit the ability to measure 7. 

3. GRAVITATIONAL-RADIATION-DRIVEN MASS TRANSFER 

We now review the basic scenario for mass transfer 
driven by angular momentum losses due to gravitational 
radiation. The gravitational-radiation-drivcn mass trans- 
fer scenario was first developed (Kraft et al. 1962; Faulkner 
1971) to explain the properties of white dwarf binaries with 
short periods, however the physics also applies to neutron 
star binaries. In this section we will give a short review 
of the derivation of the equations including the effects of 
both mass loss and the neutron star's binding energy. The 
derivation of these equations excluding the binding energy 
effect appears in many sources, such as Rappaport et al. 
(1982) and Verbunt (1993). We will make use of the fol- 
lowing notation: Mi is the neutron star's gravitational 
mass, while M2 and R2 are the mass and radius of the 
companion star. The mass ratio is denoted q = M2/M1. 
For the ultra-compact systems which are most likely to be 
of interest, the mass ratio is probably very small. 

If there is no mass transfer between the stars, the or- 
bital separation will slowly shrink due to the emission of 
gravitational radiation. The binary system loses angular 
momentum at the rate 



GB. 



-^.(1 + 9)-/^ 
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p (19) 



where P is the orbital period. As the stars move closer to- 
gether the Roche lobe radius of the secondary star shrinks 
and eventually meets the surface of the secondary, at which 
point mass will begin to be transfered to the neutron star. 
Once mass transfer begins, the size of the secondary star 
is assumed to be the same as the Roche lobe (Rl)- In 
the Paczyhski (1971) approximation, the size of the Roche 
lobe (and R2) is given by 

1/3 



R2 = Rl = 0.46a 



q 



(20) 



A + q. 

where a is the distance between the stars. The time deriva- 
tive of R2 is 

i?2 a 1 q 1 
R^^a^ 3q{l + q)' 

The secondary loses mass at the rate M2- Since the sec- 
ondary is assumed to be a non-relativistic star (such as 



(21) 




a white or brown dwarf), we do not need to make a dis- 
tinction between the gravitational and baryon masses of 
the secondary. As the secondary loses mass, its radius is 
assumed to change as 



(22) 



where Uad describes the purely adiabatic changes to the 

star that depend only on the mass loss, while {R2/R2] 

V Jth 

corresponds to changes in the star's structure due to ther- 
mal adjustments. This last term is of order \/tkh, where 
tkh is the Kelvin-Helmholtz timescalc. While in most 
cases the thermal timescale can"t be neglected, in the case 
of a cool white dwarf the Kelvin-Helmholtz timescale will 
be long enough that only terms depending on the mass 
loss rate need to be included in equation (22). 

In the most naive model of a degenerate star, = 
— 1/3. However for very low mass white dwarfs with finite 
temperature and an equation of state including Coulomb 
repulsion, the values of Uad range from -1/3 to (Deloye & 
Bildsten 2003), depending on the mass, composition and 
temperature. If the mass and composition are known, the 
uncertainty in temperature introduces an uncertainty in 
Uad of about 10%. 

For now, we will allow the most general mass transfer 
and assume that only a fraction (3 of the mass lost from 
the secondary is captured by the neutron star. As a result, 
the baryon mass of the neutron star increases at the rate 

Mo,i = -I3M2- (23) 
However, the gravitational mass of the neutron star in- 
creases at a slower rate, given by 

Ml =-(l-7)/3M2, (24) 
where the parameter 7 (discussed in the previous section) 
vanishes for newtonian stars and is approximately equal 
to the neutron star's compactness. 

Since we are allowing for non-conservative mass trans- 
fer, angular momentum will be lost. Following the work of 
Rappaport et al. (1982), we parametrize this angidar mo- 
mentum loss due to non-conservative mass transfer with 
the parameter a, so that 



27ra^ 
a{l- f3)—^M2. 



(25) 



■M p 

The parameter a encompasses the uncertainty in where 
the mass is lost. If the mass loss occurs close to the centre 
of mass (say near the neutron star) then a will be very 
small (of order g^, see Podsiadlowski, Rappaport, & Pfahl 
(2002)), while if the mass loss occurs near the low mass 
secondary, a will be close to unity. 

The liberation of energy during accretion due to the ef- 
fects discussed in this paper also introduce a mechanism 
for angular momentum loss. During the accretion process, 
the number of baryons and angular momentum transfered 
to the star are conserved, while energy is liberated at the 
rate E = (3^M2C^ in the rest frame of the star. How- 
ever, the star is orbiting the centre of mass at a distance 
oi ai = qa/ [l + q). Introducing a parameter 5, of order 
unity, the angular momentum loss due to the orbital mo- 
tion of the neutron star as the energy is liberated is 



5 



In the systems of interest, the mass ratio q is likely to be 
vc^ry small, so this term will not be of much importance. 

Another mechanism for angular momentum loss is mag- 
netic braking, where the donor star's stellar wind causes 
the donor to lose angular momentum through a coupling 
with the star's magnetic field. If the binary system is 
tidally locked to the donor's spin, this causes the binary 
system to lose orbital angular momentum (Verbunt & 
Zwaan 1981). We denote this mechanism through a term 
Jb which is independent of the mass loss rate. Estimates 
of the magnetic braking law depend on knowledge of main 
sequence stellar winds and rotation rates. This type of 
mechanism is most likely to be of importance if the donor 
was originally a main sequence star. The effects of different 
parameterizations of Jb have been explored by Rappaport, 
Verbunt, & Joss (1983). 

Since the mass accreted onto the neutron star adds an- 
gular momentum to the star, this too creates an effective 
loss of orbital angular momentum. However, since this loss 
is proportional to {Ri/af, (where R\ is the neutron star's 
radius) it is negligible. 

Combining together all the different mechanisms for loss 
of orbital angular momentum, the total angular momen- 
tum loss rate is 




J 
J 

GR \ / M \ /7 \ /B 

Since the binary's orbital angular momentum is 



J = M1M2 



Ga 



M1 + M2 



1/2 



(27) 



(28) 



the change in orbital angular momentum due to rearrange- 
ments of the mass is 



j\ _la Ml M2 1 (Ml + M2) 



J 2 a Ml M2 2 (Ml M2) ■ 



(29) 



Equations (21), (22), (24), (27) and (29) are a set of equa- 
tions that allows the mass accretion rate to be written in 
terms of the angular momentum loss rate due to various 
processes, 



M2 
M2 



1 



2A{a,(3,^,nad,q) 



GR 





where the function A{a, (3, 7, riod, q) is defined by 

-g(l-7)/3-2Y^(l-/3 + /37) 



4a(l-/3)(l + ,)-^^/37^, 



(31) 



and it should be remembered that the adiabatic index n^d 
depends on M2 = qMi. When 7 = 0, this equation re- 
duces to the equation derived by Rappaport et al. (1982) 
for non-conservative mass transfer. In the binary systems 
of interest q is small, so the term proportional to 5 is ex- 
plicitly very small. It should also be remembered that the 
constant a is proportional to q^ if the mass loss occurs near 
to the neutron star, which is the most likely situation. 

If we were to ignore changes in the secondary due to 
thermal adjustments and any angular momentum losses 
due to magnetic braking, a formula for the orbital period 



derivative can be derived. (We will the discuss the appli- 
cability of these assumptions in the next section.) Making 
use of Kepler's law, we find the following equation for the 
orbital period derivative 

967r g(l + 9)-i/3 (2t:GM^\^"' ^^^^ 



(33) 



5 B{a,f3,j,nad,q) \ c^P 
where the function B{a, (3, 7, riod, q) is defined by 

2 



B 



-A. 



1 - 3nad 

Making use of the mass function /, 

{qsmif 

^-^^(TT^' 

the orbital period derivative is 
■ Q&Ti ( f{l + q)V'^ 1 /27rGMo^^/^ 



5 V ^0 



Bsmi \ c^P 



Ml 

Mq 



(34) 



4/3 



(35) 



The predicted values of P can be written in the form 

B sm z V McT. 



(36) 



where the observationally determined constant V is de- 
fined by 



V = 18.40 X 



fl/3 

/-6 

n5/3' 



(37) 



and //Mq /_6 X 10"^ and P = P3 x 10^ s. Values for the 
constant V for the ultra-compact binaries with measured 
mass functions are given in Table 1. Since the contribu- 
tion to the period derivative from the relativistic term 7 
is quite small, we will typically require 3 significant fig- 
ures, hence the displayed accuracy in equation (37) and in 
Table 1. Although the binary mass ratio q is quite small 
for the systems of interest we have explicitly kept the fac- 
tor of (1 + q)^^^ in the equation for the period derivative 
since a precise result is required. Present observations of 
the orbital period are not precise enough to detect such a 
small period derivative in the orbit of any of the systems 
listed in Table 1, but observations over a long period of 
time may allow the detection of the period derivative. 

A similar formula for the rate that mass is lost from the 
companion can be written in the form 



A{sm.iY \MqJ yr 
where the constant is defined by 



jr= 58.1 X 



„2/3 
/-6 



P 



8/3 ' 



(38) 



(39) 



Values for the constant are given in Table 1. It has 
already been shown by Chakrabarty & Morgan (1998); 
Bildsten & Chakrabarty (2001) that the mass transfer rate 
for SAX J1808. 4-3658 is consistent with being driven by 
gravitational radiation. However, the determination of M2 
is probably no better than order of magnitude, so these 
measurements can't be thought of as providing an extra 
constraint on 7. 
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4. ASTROPHYSICAL CONSIDERATIONS 

The expression for the period derivative derived 
in the previous section depends on the parameters 
Ml , sin I, Had, a, P, S and 7. (The mass ratio q is not inde- 
pendent, since it is determined once / is measured and Mi 
and sini are specified.) This expression was derived by as- 
suming that the secondary's thermal timescale is very long, 
and that magnetic braking is not operating. We now con- 
sider some of the models for the evolution of ultra-compact 
binaries in order to cvahiate in which models there might 
be a chance of detecting 7. 

There are two main scenarios for producing a neutron 
star X-ray binary with a very short orbital period (see 
Nelson & Rappaport (2003) for a discussion). In the first 
scenario, the companion star is a main sequence star and 
will be denoted MS-NS. The alternative picture is of a 
white dwarf companion and will be denoted WD-NS. 

In the MS-NS scenario, mass transfer from a low mass 
star occurs near the time when Hydrogen is depleted in its 
core (Nelson & Rappaport 2003). Both gravitational radi- 
ation and magnetic braking (Ergma & Antipova 1999) re- 
move orbital angular momentum from the system, and the 
orbital periods can evolve down to very short timcscalcs 
(Podsiadlowski, Rappaport, & Pfahl 2002). In this type 
of system, magnetic braking could easily be as important 
as gravitational radiation in the removal of angular mo- 
mentum. Since the parameters describing magnetic brak- 
ing are not accurately known, it would not be possible to 
measure the effect of the 7 term in this type of binary 
system. 

In the WD-NS scenario the white dwarf fills its Roche 
lobe when the orbital period is very small (3 minutes for 
example) and gravitational radiation drives mass transfer 
and increases the orbital period to the range of 40 min- 
utes, reducing the white dwarf mass to very low values of 
order O.OIMq Rasio, Pfahl, & Rappaport (2000). Mag- 
netic braking is unlikely to be operating in the WD-NS 
picture, since it is most effective if the donor star has a 
large radius, and it typically assumes a convective enve- 
lope that drive a magnetic dynamo. A more difficult ques- 
tion is whether the white dwarf is being heated, through 
either being irradiated by the X-ray flux from the neutron 
star or through tidal effects. If the heating timescale is 
short compared to the mass-loss timescale, the {R2/R2)th 
term can't be discarded. Deloye & Bildsten (2003) have 
shown that observations of M2 and P for WD-NS systems 
can constrain possible combinations of WD composition 
and temperature and have demonstrated that it is possi- 
ble to use these observations to show that in some cases 
an adiabatic evolution is consistent with observations. 

If a measurement of the orbital period derivative of a 
WD-NS binary were to be made, and if one assumed that 
the evolution is adiabatic, equation (35) could be used to 
constrain a combination of Afi, sin ?', Uad, Q, ^ and 7. In 
the compact systems most likely of interest, q is proba- 
bly quite small, of order 0.01, so terms of order q^ in the 
function A (such as the term proportional to 5) can be ne- 
glected. If any mass loss take place near the neutron star 
(say through a propeller mechanism) then the parameter 
a ^ and the a term can be neglected. In the function 
A, 7 always appears in the combination g/3(l — 7), so it is 
essentially this term that we would like to measure. Since 



7 could realistically be expected to have a value near 0.2 

if we were pessimistic and supposed that (3 could be as 
small as 0.5, the term q(3{l — 7) > 4 x 10^^. We would 
then want to know Uad to an accuracy of at least ±0.001 
in order to have a hope of measuring the combination of 
/3(1 — 7), which would give an upper limit on the value of 
7- 

If the WD temperature and composition were known 
then Uad is given by Deloye & Bildsten (2003) as a function 
of Ml and sini. For a given WD temperature and compo- 
sition, a period derivative measurement (along with all of 
the other assumptions listed in this section) would yield a 
constraint on a combination of Mi, sinz, and /3(1— 7). Fur- 
ther constraints on the inclination angle would strengthen 
this constraint. For example recent observations of SAX 
J1808.4-3658 by Wang et al. (2001) have made use of a 
reddening distance technique and models of X-ray heated 
accretion disks to constrain the binary inclination, depend- 
ing on the assumed distance to the binary system. 

In the equations for the period derivative for mass trans- 
fer induced only by gravitational radiation with adiabatic 
changes to the white dwarf, the value of the Uad plays the 
most important role (along with the inclination angle and 
the neutron star mass) in determining the binary's evo- 
lution. Ideally, observations of such as system ignoring 
the relativistic 7 term could be used to constrain ranges 
of values of nad consistent with an adiabatic evolution for 
various possible values of sin i and Mi . Once these values 
are determined, and expanding equation (31) to lowest 
order in q (assuming that the mass loss occurs near the 
neutron star) one finds the limit on 7, 
^3 l/(5 + 3n„d) 

(l-3n„,)PxlO-i^ (l + g)V^ ( M^X"\,^. 
2 p smi \Mq) f ' 

where we have used the fact that /? < 1. Although a num- 
ber of assumptions have been made in deriving this limit, 
inequality (40) does have the potential to provide an in- 
teresting limit on 7. Since the mass ratio q is expected to 
be small (of order 0.01), a meaningful limit will only occur 
if the terms proportional to 1/ q are similar in magnitude, 
so that their difference is of order q. Since V is of order 
unity, only orbital period derivatives of order 10~^''s/s will 
provide an interesting limit on 7, given all of our assump- 
tions. From equation (40) it can be seen that if P is too 
small, the maximum allowed value of 7 will be so large as 
to be an uninteresting limit. If however, P is too large, 
the maximum allowed value of 7 will be negative, implying 
that our assumptions about the physics of the binary are 
invalid. 

As an example, consider the situation where the com- 
panion star is a Helium white dwarf with mass O.O2M0 
and an adiabatic index of nad = —0.22, consistent with a 
temperature near IQ^ K as calculated in the models of De- 
loye & Bildsten (2003). With data for the mass function 
for XTE J1751-305, inclination angles with sini > 0.7 give 
neutron star masses larger than 1.4 Mq. Given this range 
of inclination angles, it is easy to use equation (40) to find 
the maximum allowed value of 7 for any observed P. These 
values are shown in Figure 3. This figure can be inter- 
preted by noting that if the binary's inclination angle were 
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Fig. 3. — Plots of the maximum allowed values of 7 versus binary inclination angle for XTE J1751-305, assuming a Helium white dwarf 
companion with M2 = 0.02Mq and an adiabatic index of Uad = —0.22. In this plot, each curve corresponds to a value of the orbital period 
derivative in units of 10~^*s/s. 



found to be sin i = 0.9, then a value of P = 3.9 x IQ-'^^s/s 
would be consistent with our assumptions, while a larger 
period derivative (say 4.1 x 10"-^^) would yield a negative 
7 and is inconsistent with our assumptions. 

5. CONCLUSIONS 

In this paper we have discussed the effect on the orbital 
period of changing a neutron star's binding energy through 
accretion. We have shown that it is possible, in principle, 
for an observation of an orbital period derivative and with 
some knowledge of the donor star to place constraints on 
possible values of the neutron star's mass and the binding 
energy parameter that depend on the binary's inclination 
angle. However, there are a number of uncertainties in- 
herent in our method. We will now discuss these potential 
problems. 

Our calculation assumes that the binary's orbit is given 
by newtonian physics and have not considered post- 
ncwtonian (PN) corrections. However, the 1st PN cor- 
rections enter at the level of v'^/c^ where v is the orbital 
velocity. Since we are considering wide (from a relativist's 
point of view) binaries with separations of the order of 
10'' km, the 1st PN corrections will be about 1000 times 
smaller than the orbital corrections due to the change in 
binding energy discussed in this paper. 

The model of mass transfer considered in this paper 
makes use of the Roche model, ignores the spin angular 
momentum of the secondary, and calculates the angular 
momentum loss by gravitational radiation in the point par- 
ticle approximation. The validity of these assumptions has 
been examined by Rezzolla et al. (2001) for the case of a 
white dwarf primary. Rezzolla et al. (2001) do an accu- 
rate numerical integration to compute the structure of a 
semi-detached compressible co-rotating binary in hydro- 
static equilibrium. They find that Paczynski's approxima- 
tion for the Roche lobe radius and the expression for the 
total angular momentum are correct to about 1% if the 
secondary star's spin angular momentum is included. In 
their computations, they also compute the rate that an- 
gular momentum is lost when the secondary's structure 
is taken into account and find that the point-particle ap- 
proximation introduces an error of about 1%. While these 
corrections are quite small, the relativistic corrections due 
to the neutron star's structure are also quite small. For 
this reason it would be desirable to consider more realistic 
models such as those computed by Rezzolla et al. (2001). 

Larger uncertainties lie in the assumed physics of the bi- 
nary system. We need to restrict our attention to systems 
where the donor star is a white dwarf, and the angular mo- 
mentum loss is strictly through gravitational radiation. It 



is also important that the donor star is not rapidly heated 
so that it is possible to assume that the donor's radius 
responds adiabatically to mass loss. 

If any mass is lost from the binary system, the angular 
momentum lost could potentially drive a period deriva- 
tive correction which is larger than the relativistic effect 
considered here. If the mass loss occurs near the neutron 
star, (through the propeller mechanism, for example) the 
effect is not large enough to be of importance for very low 
mass companions. If mass is lost from a region near the 
companion, then the large angular momentum loss would 
make it impossible to measure the relativistic effect. 

The gravitational-radiation-driven angular momentum 
loss discussed in this paper is a continuous process, while 
the mass transfer onto the star is highly episodic in that 
matter can be stored in the accretion disk for years be- 
fore being transfered onto the neutron star. Since our 
treatment assumes a continuous transfer of mass, mea- 
stiremcnts of the orbital period derivative would have to 
be measured over the course of many years before there 
could be any hope of measuring the relativistic parameter 
7- 

The orbital frequencies typical for the binaries of inter- 
est lie within the bandwidth of the Laser Interferometer 
Space Antenna (LISA). The dimensionless gravitational- 
wave amplitude h (Thorne 1998) measured by LISA can 
be written as 

where D is the distance to the binary. Using the measured 
binary parameters for SAX J1808. 4-3658, we find that the 
gravitational- wave amplitude is h < 10~^^, which is well 
below the "confusion noise" produced by the population 
of galactic white dwarfs (Nelemans et al. 2001). The low 
amplitude of the signal is mainly due to the small mass 
ratio, so we expect that the binaries of the type discussed 
in this paper will not be detectable by LISA, unless new 
data analysis techniques are developed which can take ad- 
vantage of the known positions of the binaries. 

In order for the effects of a changing relativistic mass 
deficit to be observable, very small changes in the binary's 
orbital period must be measured. In addition, the rate of 
mass transfer must be known in order to calculate the mag- 
nitude of the relativistic effect. For this reason, we have 
examined the scenario of mass transfer from a fully de- 
generate dwarf driven by gravitational radiation, since the 
mass transfer rate is related quite simply to the binary's 
parameters. We have focused our attention on the ms 
spin-period neutron stars since the prospects for measur- 
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ing the orbital period for these systems seems most promis- 
ing. However, the orbital effect of a changing mass deficit 
also holds for non-rotating stars. In fact, the effects are 
much simpler in the case of a non-rotating star, since the 
effects of torquing due to mass accretion are completely 
negligible. 
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Table 1 

Observational data for accreting ms pulsar systems and values of the constants O and J^. 



Binary System 


Us (Hz)- 


i'.(Hzs-i) 


P (s)b 


/ {MqT 






Reference 


SAX JI808.4-3658 


401 


< 7 X 10-13 


7249.119 


3.7789 X 10-5 


2.274 


3.3 


Chakrabarty & Morgan (1998) 


XTE J 175 1-305 


435 


< 3 X 10"" 


2545.3414 


1.2797 X lO"'^ 


4.210 


5.7 


Markwardt et al. (2002) 


XTE J0929-314 


185 


-9.2 X 10"" 


2614.746 


2.7 X 10"^ 


1.1 


0.4 


Galloway et al. (2002) 


XTE J1807-294 


191 




2404.2 








Markwardt et al (2003) 


XTE J1814-338 


314 




15390 








Strohmayer et al (2003) 



-Pulsar spin frequency 

^Orbital Period 
■^Mass function 

'^V is defined in equation (37). 
'^J^ is defined in equation (39). 



